Direct Sampling of Negative Quasiprobabilities of a Squeezed State 



T. Kiesel and W. Vogel 

Arbeitsgruppe Quantenoptik, Institut fiir Physik, Universitlit Rostock, D-18051 Rostock, Germany 

B. Hage 

ARC Centre of Excellence for Quantum-Atom Optics, Department of Quantum Science, 
The Australian National University, Canberra, Australian Capital Territory 0200, Australia 

R. Schnabel 

Institut fur Gravitationsphysik, 
Leibniz Universitdt Hannover and Max-Planck-Institut fiir Gravitationsphysik (Albert-Einstein-Institute), 
Callinstrasse 38, 30167 Hannover, Germany 

Although squeezed states are nonclassical states, so far, their nonclassicahty could not be demonstrated by 
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of so-called nonclassicahty quasiprobabilities. The negativities of these quantities turn out to be necessary 
and sufficient for the nonclassicahty of an arbitrary quantum state and are therefore suitable for a direct and 
general test of nonclassicahty. We apply the method to a squeezed vacuum state of light that was generated by 
parametric down-conversion in a second-order nonlinear crystal. 
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Introduction. In quantum optics and quantum informa- 
tion science, the notion of nonclassicality describes the dis- 
tinguished difference between classical and quantum physics. 
Here, a quantum state is referred to as nonclassical, if one is 
not able to model the outcomes of experimentally measured 
optical field correlation functions by classical electrodynam- 
ics. Considering solely pure states, the famous coherent states 
I a) are the only classical states according to this notation, 
which makes them the closest analogue to the classical os- 
cillator [1]. Sudarshan [2] and Glauber [3| showed that the 
density operator of an arbitrary quantum state can formally be 
written as a statistical mixture of coherent states, 

J (faP{a)\a){a\. (1) 

If P{a) resembles the properties of a classical probability 
function, the state is simply a classical mixture of the (clas- 
sical) coherent states, e.g. a thermal state. In general, the P 
function may attain negative values - often in connection with 
a strongly singular behavior. In such cases the corresponding 
quantum state is referred to as a nonclassical one (lU. 

The main problem of this definition of nonclassicality lies 
in the singularities of the P function, which definitely prevent 
the experimental reconstruction of P{ct). Only for special 
quantum states one may approximately obtain this quasiprob- 
ability 0. Therefore, different criteria for the detection 
of nonclassicality have been developed. Some of them are 
simple, such as squeezing f6|, classical limits on probabili- 
ties |7| or negativities in the Wigner function |8|, but they 
are only sufficient for nonclassicality. Others are necessary 
and sufficient, but they consist of an infinite hierarchy of in- 
equalities EllTOl. Recently, nonclassicality quasiprobabili- 
ties have been introduced, which provide a complete and sim- 
ple method for the verification of nonclassicality tlli : For 



any nonclassical state, there exists a regular nonclassicality 
quasiprobability, whose negativities unambiguously reflect its 
nonclassicality. In flZl, the experimental applicability, as a 
matter of principle, was demonstrate on a nonclassical but less 
problematic state, which also had a negative Wigner function. 

In this Letter, we prove the nonclassicality of a Gaussian 
squeezed state by reconstructing negative quasiprobabilities 
from data taken by a balanced homodyne detector. Squeezed 
states are of special interest in this context, since their com- 
monly used quasiprobabilities, such as the Wigner function, 
are nonnegative and still satisfy the properties of classical 
probability densities. We avoid any Fourier transformation 
of the data, which was used in [12|, and present a method of 
direct sampling of nonclassicality quasiprobabilities from the 
measured data. For this purpose, we use the concept of pattern 
functions ifTSll . which provide an estimate of the quasiproba- 
bility together with its variance. This method applies to the 
experimental characterization of nonclassicality of any quan- 
tum states, the only limitation being statistical uncertainties. 

Quasiprobabilities of squeezed states. Squeezed states are 
prominent examples of nonclassical states, which can be eas- 
ily generated in the laboratory. Although nonclassicality is 
defined by negativities of the P function, its general verifi- 
cation by negativities of any commonly used quasiprobability 
distribution is impossible. For instance, the Wigner function 
of a squeezed state with quadrature variances Vx and Vp reads 
as 

clearly being a Gaussian. In contrast, the P function may 
formally be written as 

P^^{a) = e {l^ + in^ vx-vp f°a^;,S(a), (3) 
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representing one of the most singular representations of a 
quantum state, with infinitely high orders of derivatives of the 
(5-distribution. Moreover, the s-parameterized quasiprobabili- 
ties 1 14 1 of a squeezed state are either Gaussian (and nonneg- 
ative) or strongly singular Based on a simple condition for 
the characteristic function of the P function |fT5|. the nonclas- 
sicality can be easily verified |16|. However, this condition 
is sufficient only and the question remains if there exists any 
well-behaved quasiprobability, allowing a complete character- 
ization of the nonclassicality of squeezed states by the failure 
of being interpreted as a classical probability. 

Nonclassicality quasiprobabilities. The starting point of 
our discussion is the characteristic function of the P function, 

$(^) = (: el^^^-r^ :) = (^(,m3:[<^^K{P)-^ /■A-^f.m^ /■i ^ (4) 

with x{ip) being the quadrature operator of the optical field at 
phase ip. In order to obtain a regular phase-space distribution, 
we filter in the form 



$n(/3) = <f(/3)f^»(/3). 



(5) 



The filter Vl^^ji) has to satisfy the following conditions to be 
useful for nonclassicality detection 111]: 

1. The filtered characteristic function <i>o(/3) should be in- 
tegrable for an arbitrary quantum state $(/3), such that 
its Fourier transform - the nonclassicality quasiproba- 
bility - exists as a regular function. 

2. Negativities in the Fourier transform shall unambigu- 
ously be due to the nonclassicality of the state. Con- 
versely, the nonclassicality quasiprobability shall be 
nonnegative for any classical state. This requkes that 
the filter Sl^ has a nonnegative Fourier transform. 

3. If the width parameter approaches infinity, the filtered 
characteristic function <I'o(/3) should converge to the 
characteristic function of the P function, $(/?). This 
can be realized by the conditions ^w{P) = ^liP/w) 
andr2i(0) = 1. 

4. The filter should be nonzero everywhere, fi^ (/?) ^ 0, 
such that no information about the quantum state is lost 
due to the filtering in Eq. (|5]l. 

Under these conditions, the nonclassicality quasiprobability 
is defined as the Fourier transform of the filtered characteristic 
function. 
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ctl3* — a* I 



(6) 



Negativities of the quasiprobability are signatures of nonclas- 
sicality of the state and the negativity of its P function, since 
the filter is constructed in such a way that it does not introduce 
negativities by itself. In the present work, we construct a filter 
from the autocorrelation of the function uj{/3) = exp(— 1/3|"'), 



1 



(7) 



the normalization constant is chosen to obey r2i(0) = 1. 
The width is introduced via flwiP) — fli{f3/w). This filter 
satisfies all criteria mentioned above, for the proof see fTT]. 

Derivation of a pattern function. Pattern functions pro- 
vide an efficient technique to directly estimate a physical 
quantity together with its uncertainty. From balanced homo- 
dyne detection, we obtain quadrature values Xj{Lp) measured 
for certain phases Lp. They obey the quadrature distributions 
p{x; (p), which satisfy J p{x; ip)dx = 1. The quadratures are 
normalized such that the vacuum quadratures have a variance 
Vvac = 1- Now the nonclassicality quasiprobability Pn{ct) 
with a certain width parameter w is written as the statistical 
mean of the pattern function /o {x, (p; a, w), averaged over the 
quadrature distributions: 

p(x;ip) 



dx dip' 
Jo 



-fn{x,(p;a,w). (8) 



For this purpose, we note that due to Eq. Q, the characteristic 
function of the P function of the state can be calculated from 
the quadrature distribution via 



dx p (x; arg(/3) - f ) e'l^l^gl^l'/^. (9) 



It is convenient to rewrite the integral in Eq. (j6]l in polar coor- 
dinates /3 = 5e"''. Here, we restrict <pto [0, tt) and extend b to 
(—00,00). Then we obtain 



db / d(p|6|e2*l°l''"'"("s(«)-¥')$(-^e*¥') 



xa^ibe'"^). (10) 



The filter is chosen to be independent of the phase, 
i.e. r2uj(6e"'') = fliu{b). Now we insert Eq. (|9| and obtain 



dx / dip 
00 Jo 



p{x;lp) 



^g2i|a|fcsin(arg(a)-<^-2)l7^(6). 

This equation defines the pattern function 



(11) 



fn{x,(p;a,w) 



Jfl 1^1 ^ibx ^2i\a\bsin{&rg(a)-ip-^) 



(12) 



which has to be used in Eq. ([8]). 

Equation ^ gives rise to the following interpretation: Sup- 
pose we have measured N quadrature-phase pairs {xi,(pi), 
whose joint probability distribution is ^p{x;(p). Here 
the phases (p are assumed to be uniformly distributed in 
[0, tt), while the quadratures obey the quadrature distribution 
p{x; ip), conditioned on the value of the phase pj. Then the 
quasiprobability Pn{a) can be calculated as the expectation 
value of the pattern function /a(a;, ip; a, w). For experimental 
data, we replace the expectation value by its empirical esti- 
mate. 



Pn{a) 



1 ^ 
N ^ 

i=l 



/o(a;j,i^i;a,w;). 



(13) 
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Its variance can be obtained naturally as the mean square de- 
viation of the numbers fn{xi, ipi; a, w). 

If the phases, at which quadratures are measured, are 
scanned in [0, tt] or drawn randomly from a uniform distri- 
bution, one can calculate the nonclassicality quasiprobability 
directly as the empirical mean of the pattern function. This 
mean is taken over all pairs {xi, ipi) of quadrature and phase. 
In our experiment, we only obtained quadratures at 21 fixed 
phase angles. In this case, one may not simply replace the 
integral over the phase (/? by a sum: This leads to systematic 
deviations, since the pattern function is varying rapidly with 
respect to the phase, in particular if \a\ becomes large. For a 
detailed discussion and solution of this problem, we refer to 
the supplemental material I J 71 . 

Experimental set-up. The squeezed vacuum states of 
light were generated by type-I degenerate parametric down- 
conversion (optical parametric amplification, OPA) inside 
an optical resonator. The latter was a standing wave res- 
onator with a line width of 25 MHz containing of a non- 
critically phase matched second-order nonlinear crystal (7% 
Mg:LiNb03). The OPA process was continuously pumped by 
50 mW of second harmonic light yielding a classical power 
amplification factor of six. Both, the length (resonance fre- 
quency) of the resonator as well as the orientation of the 
squeezing ellipse were stably controlled by electronic servo 
loops. With this setup we directly measured a squeezed vari- 
ance of -4.5 dB and an anti-squeezed variance of H-7.2 dB with 
respect to the unity vacuum variance. From these measure- 
ments we inferred an overall efficiency of 75% and an initial 
squeezing factor of -8.2 dB. 



OPA 




Figure 1. Simplified sketch of the experimental setup. A spatially 
filtered continuous-wave field at 1064 nm served as a local oscillator 
(LO) for balanced homodyne detection (BHD) and a phase-locked 
second harmonic field at 532 nm as the pump for the parametric 
squeezed light source (OPA). MC: spatial mode cleaner, PZT: piezo- 
electrically actuated mirror for adjusting the quadrature amplitude 
phase of the BHD. 

The quadrature amplitudes of the squeezed state were mea- 
sured by balanced homodyne detection (BHD). The visibility 
of the squeezed field and the spatially filtered (MC, Fig. [T]i 
local oscillator was 98.9% and was limited by OPA crystal 
inhomogeneities. The quadrature phase of the BHD was ad- 
justed by servo loop controlled micro-positioning of steering 
mirrors in one of the optical input paths. The photo-electric 
signals of the two individual BHD-photodiodes were electron- 
ically mixed down at 7 MHz and low pass filtered with a band- 



70 r 

60 

_ 50 

w_ 

40' 

C 

o 

1 30 
20 



10 



1.2 



1.4 

w 



1.6 



1.8 



Figure 2. Absolute value of the significance E(ui) of the negativity 
of the quasiprobability versus the filter width w. 



width of 400 kHz to address a mode showing good squeezing 
and a high detector dark noise clearance of the order of 20 dB. 
The resulting signals were fed into a PC based data acquisition 
system and sampled with one million samples per second and 
14 bit resolution and finally subtracted yielding the quadra- 
ture amplitude data. For a more detailed description of the 
main parts of the setup we refer to ifTSl [T9]| . 

Experimental results. The examined squeezed vacuum 
state is characterized by the variances 14 = 0.36 and Vp ~ 
5.28. We acquired IQ^ quadrature values for each of the 21 
quadrature phases, the latter being equally spaced in [0,7r]. 
The values of the quasiprobability Pn [a) as well as their stan- 
dard deviation (T(Po(a)) are estimated from the pattern func- 
tion as given in Eq. ( [T5] l. The filter width is chosen such that 
the significance of the negativity is optimized. Our figure of 
merit is defined as 



Pn{a) 



(14) 



with Tj{w) being negative if is negative for some a. 

By construction of the quasiprobability, a significant negativ- 
ity clearly indicates nonclassicality of the state. The larger the 
absolute value of Yj{w), the larger is the significance of the 
negativity. This quantity can be optimized with respect to w. 
In Fig|2] we show the dependence of the significance on the 
filter width w. The larger the filter width, the larger the non- 
classical effects of the state are manifested in negativities, but 
the larger also the statistical uncertainty grows. Therefore, an 
optimum width exists, which is achieved here for w = 1.3. 

Figure |3] shows the nonclassicality quasiprobability ob- 
tained from the experimental data. We find that along the axis 
of Im(Q:), the quasiprobability oscillates and becomes clearly 
negative. This uncovers the nonclassicality of the squeezed 
state in a general sense, beyond the phase-sensitive reduction 
of the quadrature variance. It also includes the information 
on other kinds of effects, such as higher-order squeezing of 
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Figure 3. Nonclassicality quasiprobability for a squeezed vacuum 
state. The data is directly sampled from our balanced homodyne 
data and clearly shows negative values. 



0.6 




-2-1 1 2 

lm(a) 



Figure 4. Cross section of the nonclassicality quasiprobability for 
our squeezed vacuum state. Note, that the uncertainty in the data is 
less than the line width chosen here. 



Conclusions. In our work, we introduced a method for the 
direct sampUng of nonclassicahty quasiprobabihties of arbi- 
trary quantum states from measured quadrature ampHtudes. 
By applying our method to a squeezed state, whose P func- 
tion belongs to the most singular ones, we experimentally ver- 
ified nonclassicality in its general sense, i.e. through negative 
quasiprobabihties. Our method is not only capable of estimat- 
ing the negativity of the quasiprobability, but also its statistical 
uncertainty, in a surprisingly simple manner. 
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the types considered in ifTOll and ||20ll , and others, see |fT6ll . In 
Fig.|4] we show a cross-section of Fig. [3] along this axis. We 
clearly observe distinct negativities. The standard deviation 
is less than 1.1 x lO^'^ for all points and therefore covered 
by the width of the line. We also calculated a systematic er- 
ror due to the finite set of examined phases, being less than 
3.6 X 10"^ for all points along this axis IfTTl . attains 
the minimum at a = 0.9i with Psi(amin) = —0.05989 and 
c(-Pa(Q;min)) = 0.9 X lO^'^, therefore leading to a signifi- 
cance of |E| ~ 69 standard deviations. Hence, this is a very 
clear demonstration of the nonclassicality of the examined 
state by means of negativities of a nonclassicality quasiprob- 
ability, which is not possible for commonly used quasiproba- 
bihties such as the Wigner function. 
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Supplementary Material 

Numerical evaluation of the pattern function. The pattern 
function is given by 



1^1 ^i6z^2i|a|bsin(arg(a)-y-f ) 



xe'''/2r2„(6). 



(15) 



Since it has to be calculated for a large set of data points, its 
fast evaluation is of great relevance. First, we note that it can 
be regarded as a function of two parameters. 



d6^e*''«e'''/2r!„(6), (16) 



which is connected to frt{x, ^\ a, w) via 

fa(x, LP] a, w) = x(a;+2|a| sin(arg(Q!)-v3-7r/2); uj). (17) 

The data points only enter in x via the argument ^, while w is 
arbitrary but fixed for all data points. Therefore, it is conve- 
nient to determine the function xi^'^w) in advance. This can 
be done by Fourier techniques: Since il^(fo) decays very fast 
with h, we may approximate it by setting it to zero for all h 
with \h\ > be, with some be being sufficiently large. Note that 
be scales with the width w, since fl^^ (/?) is just a fixed function 
with a scaled argument /3 — >■ f3/w. With this assump- 
tion, is ^ bandlimited function, i.e. its Fourier trans- 
form has bounded support. Therefore, the Nyquist-Shannon 
sampling theorem holds [l]: 



x(?;w) 



oo 

E 



(18) 



Hence, the Fourier integral (16 1 has to be evaluated at a dis- 
crete set of points = Afterwards, the data points are 
inserted via ^ in Eq. (18 1. In our calculations, the support 
of flw{/3) was taken to be be = 4w, where the filter n,^{be) 
became less than 10~^^. Furthermore, we evaluated 256 co- 
efficients x{jf-',w) for a numerical accuracy of 3.5 x 10^^ of 

Sampling from a discrete set of phases. The points of 
the nonclassicality quasiprobability can be obtained from the 
quadrature distributions by 



Pnia) 



f{x, ip; a, w)dtpdx. 



(19) 



Let us now consider the case that one has measured quadrature 



distributions at a certain set of phases {pkl^Zv 

As a first idea, one might replace the integral over the phase 
(y9 by a Riemannian sum: 



k=l 



p{x;ipk)fnix,ipk;a,w)dx. (20) 



However, this leads to significant systematic deviations, such 
that does not appear to be integrable. This can be seen 
as follows: The function xi£,', w) is the Fourier transform of 
some integrable function, see Eq. (16 1. Due to the Riemann- 
Lebesque lemma |2|, when ^ approaches infinity, x(^;w) 
tends to zero. In terms of a, this can be achieved if |a| — > oo, 
but arg(a) — (pk 7r/2 for all (pk- However, what happens 
if the latter condition is not satisfied? Then, fn{x, ipk; a, w) 
appears to be independent of |q;| and does not approach zero 
for large a, and the same holds for Psi(a)- Therefore, when 



one reconstructs Psi(a) via Eq. (20 1, for all arguments of a. 



which match one of the phases examined in experiment, the 
result does not approach zero but tends to a finite value. 
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Figure 5. Application of the reconstruction formulae ^19\ (solid line) 
and l |20^ (dashed line) for a vacuum state. It is clearly seen that the 
latter graph does not tend to zero for large a. The inset shows the 
dependence of Pn(/3) on the phase, with |a| — 10. 

This effect can already be observed for phase-independent 
states, such as vacuum. Figure [5] shows the exact curve ( 19 1 
(solid line) and the approximation ( [20| (dashed line) for a vac- 
uum state. We assume that the state has been examined at 21 
phases, and plot the cross section of Pq along one of them. 
It is clearly seen that the approximation ( pO] ) does not tend to 
zero for large a. The inset shows the dependence on the phase 
of a, for |a| ~ 10. We observe an oscillating function with 
21 maxima, which do not vanish for larger a. This artefact 
is caused by the fact that the pattern function fn{x, (p; a, w) 
strongly oscillates with (p when a increases. Therefore, the 
phase integral in Eq. ^T9\ cannot be properly approximated by 
a sum. 

In the following, we assume that only the quadrature dis- 
tributions p{x; (p) do not change significantly within the inter- 



vals [ipk 



2N, 



2N, 



] for all fc = 1, 



, N^. More il- 



lustratively, we assume that our set of quadrature distributions 
measured at the phases ipk suffices to have complete infor- 
mation about the quadrature distributions at phases ip which 
have not been measured. If this were not the case, one had 
to increase the phase resolution in order to obtain more in- 
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Figure 6. Systematic error due to a discrete set of examined phases. 



following way 

Pn{a) - 



E 

fc=i 



OO nTT 



f{x, ip; a, w)dipdx 



p[x; ipk + v) 



/(x, Lpk + w)d(pdx 



2N^ 



k=l 



IT 



, ^ I I J f{x,(pk + (p;a,w)d(pdx 



2N^ 



1 f 

= P{x;ipk)f{x,(pk;a,w)dx, (21) 

k—1 

where the modified pattern function is defined as 



f{x,(pk;a,w) 



2N^ 



2N^ 



f{x,tpk + (p;a,w)d(p. (22) 



The difference between the first and the last line in Eq. (21 
is a systematic error, caused by measuring at a discrete set 
of phases. Figure [6] shows this quantity for a squeezed vac- 
uum state with the measured variances = 0.36 and Vp — 
5.28 along the axis where negativities of the nonclassicality 
quasiprobability are found. There, the systematic deviations 
are approximately one order of magnitude below the statis- 
tical standard deviation, which has been given in the Letter. 



The modified pattern function (22 1 can be interpreted as the 



original pattern function, averaged over the phase ip, which 



is uniformly distributed over the interval 



2N.^ ' 2N^ 



We 



realize this by drawing a random number ip in this interval 
for each quadrature point Xj{(pk), which has been measured 
at phase ipk, and add this random number (p to the phase (pk- 
In this way, ipk represents a random phase, whose mean is 
centered at ipk, but whose values are drawn from an interval 
of the length of the difference of measured phases. With this 
Monte-Carlo-like approach, we efficiently calculated the non- 
classicality quasiprobabihties. 
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